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NONLINEAR  ANALYSIS  OF  RELATIVISTIC  HARMONIC  GENERATION  BY 

INTENSE  LASERS  IN  PLASMAS 

I.  INTRODUCTION 

The  interaction  of  ultra-intense  laser  pulses  in  plasmas  is  rich  in  a  variety  of  phenom¬ 
ena  [l]-[22].  These  phenomena  have  gained  particular  relevance  due  to  the  development 
of  compact  lasers  which  produce  ultrashort  pulses  (<  1  ps)  at  ultrahigh  powers  (>  1  TW) 
and  intensities  (>  10^*  W/cm*)  [2],  [23].  For  example,  the  production  of  harmonic  ra¬ 
diation  may  occur  by  several  mechanisms.  At  modest  intensities,  lasers  interacting  with 
neutral  gases  have  been  observed  to  produce  coherent  harmonic  radiation  at  well  past  the 
53'*  harmonic  due  to  atomic  effects  [13].  At  ultrahigh  intensities,  a  gas  is  readily  ionized 
and  the  effects  of  the  free  plasma  electrons  become  exceedingly  important.  The  ioniza¬ 
tion  process  itself  results  in  electron  currents  which  can  produce  harmonic  radiation  [14], 
[15].  In  a  fully  ionized  plasma,  harmonics  can  be  produced  by  (i)  relativistic  harmonic 
generation  in  the  forward  direction  (the  propagation  direction  of  the  incident  laser)  [1],  [8], 
[17]-[20],  (ii)  stimulated  backscattered  harmonic  generation  [1],  [21]  and  (iii)  incoherently, 
by  nonlinear  Thomson  scattering  [1],  [22].  Recent  experiments  on  these  processes  have  also 
been  perfoir>'red  [2],  [15]-[18].  This  paper  discusses  the  details  of  process  (i),  relativistic 
harmonic  generation  in  fully  ionized  plasmas.  The  growth,  dephasing  and  saturation  of 
the  harmonics  will  be  analyzed. 

An  intense  laser  field  interacting  with  a  plasma  induces  transverse  currents  associated 
with  the  quiver  motion  of  the  electrons.  For  ultrahigh  intensities  and  linear  polarizations, 
the  induced  plasma  current  becomes  highly  relativistic  and  nonlinear,  resulting  in  the 
generation  of  coherent  harmonic  radiation  in  the  forward  direction  (see  Fig.  1).  The 
transverse  quiver  momentum,  p^,  of  an  electron  in  a  ID  laser  field  is  given  by  p,  =  mocao, 
where  ao  =  eAo/wioc^  is  the  normalized  vector  potential  of  the  incident  laser  field.  The 
quiver  velocity,  v,,  is  given  by  v,  =  cao/7j.,  where  yj.  ~  (H-Oq)^^’  is  the  relativistic  factor 
associated  with  the  transverse  electron  motion.  Consider  a  linearly  polarized  incident 
laser  field  of  the  form  ao(2,t)  =  Oq  cos^oe*>  where  9q  =  koz  —  uot  and  ko  and  ti;o  are  the 
wavenumber  and  frequency  of  the  incident  laser  field,  respectively.  In  the  mildly  relativistic 
limit,  dg  "C  1,  the  quiver  velocity  is  sinusoidal.  At  ultrahigh  intensities,  dg  >  1,  the  quiver 
velocity  contains  a  spectrum  of  harmonic  components.  This  nonlinear  electron  quiver 
motion  leads  to  the  generation  of  relativistic  harmonic  radiation  [1],  [8],  [17]- [20].  The 
MMiMci^itVpproved  November  5,  1992 
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laser  strength  parameter,  do,  is  related  to  the  intensity,  /o»  of  the  incident  laser  field  by 

do  8.5  X  10-'®Ao[/im]Ji/*[W/cm^],  (1) 

and  to  the  laser  power,  Po,  by  Po[GW]  ~  21.5(doro/Ao)^»  where  ro  is  the  spot  size  of  the 
Gaussian  transverse  profile  and  Ao  =  27r/fco  is  the  incident  laser  wavelength.  For  Ao  — 
1  ^m,  ultrahigh  intensities  Jo  ^  10^*  W/cm*  imply  dg  ^  1  and,  hence,  highly  nonlinear 
and  relativistic  electron  motion.  Such  intensities  are  currently  available  from  compact 
laser  systems  based  on  the  method  of  chirped-pulse  amplification  [2],  [23]. 

Relativistic  harmonic  generation  was  first  described  and  analyzed  by  Sprangle,  et  al. 
[8].  In  Ref.  [8],  the  independent  variables  ^  =  z  —  ct  and  t  =  t  were  used  along  with  a 
quasi-static  plasma  response.  Expressions  were  derived  for  the  growth  of  the  harmonic 
radiation  in  the  linear  regime  in  which  the  harmonic  field  amplitude  is  proportional  to  the 
laser-plasma  interaction  distance,  L.  These  expressions  are  valid  for  interaction  distances 
less  than  the  phase  detuning  length,  ct^.  Use  of  the  variables  (^,t  along  with  the  quasi¬ 
static  approximation  is  not  adequate  to  accurately  describe  saturation  of  the  relativistic 
harmonics  by  phase  detuning.  More  recently,  by  using  the  variables  ^  =  z  —  and 
T  =  t,  where  is  the  phase  velocity  of  the  incident  laser  field,  along  with  the  quasi¬ 
static  plasma  response,  the  authors  accurately  analyzed  the  saturation  of  the  relativistic 
harmonics  by  phase  detuning  [Ij.  Independent  analyses  of  saturation  of  the  third  harmonic 
by  phase  detuning  have  also  been  recently  performed  [19],  [20]. 

This  paper  is  organized  as  follows.  The  remainder  of  the  introduction  describes  the 
basic  physics  of  third  harmonic  generation  using  a  simplified  ID  model  in  the  mildly 
relativistic  limit,  Oq  <$:[  1.  The  importance  of  collective,  space-charge  and  detuning  effects 
are  discussed.  It  is  shown  that  the  self-consistent  plasma  density  response  significantly 
reduces  the  source  current  driving  the  harmonic  radiation.  Phase  detuning  places  a  severe 
limit  on  the  growth  of  the  harmonic  radiation.  In  Section  II,  a  ID  nonlinear  model  valid  for 
ultrahigh  intensities,  dg  >  1,  is  formulated  and  used  to  study  the  generation  of  coherent 
radiation  at  odd  harmonics.  A  general  expression  for  the  nonlinear  index  of  refraction 
and  the  dispersion  relation  for  a  laser  field  in  the  limit  dg  >  1  is  presented.  The  collective 
plasma  response  is  included  self-consistently  and  the  saturation  amplitude  of  the  harmonics 
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by  phase  detuning  is  calculated.  This  is  done  in  the  long  pulse  limit,  ctl  »  Ap,  where 
Ti  is  the  incident  laser  pulse  duration  and  Ap  =  2iTcfup,  where  utp  is  the  electron  plasma 
frequency.  For  a  T£,  =  1  psec  laser  pulse,  cr^  >  Ap  implies  no  >  10**  cm“’,  where  no 
is  the  ambient  plasma  density.  In  Section  III,  the  effects  of  diffractive  spreading  of  the 
radiation  fields  are  determined  using  a  3D  model  in  the  mildly  relativistic  limit,  Og  'C  1. 
Third  harmonic  generation  from  semi-infinite  and  finite  slab  plasmas  are  analyzed.  The 
effects  of  transverse  gradients  in  the  initial  plasma  density  are  discussed,  which  lead  to  the 
generation  of  radiation  at  even  harmonics.  A  conclusion  is  presented  in  Section  IV. 

A.  Quiver  Model,  Collective  Effects  and  Detuning 

The  process  by  which  relativistic  effects  produce  coherent  harmonic  radiation  may  be 
understood  by  considering  a  simplified  “quiver”  model,  which  includes  only  the  effects  of 
the  relativistic  electron  quiver  motion.  Other  effects,  such  as  the  plasma  density  response, 
will  be  discussed  below.  A  linearly  polarized  laser  field  will  be  assumed,  00(2,  t)  =  do  cos  ■ 
In  the  quiver  model,  the  transverse  plasma  current  is  J,  =  —  enou,,  where  v,  is  the  relativis¬ 
tic  electron  quiver  velocity  (discussed  above)  and  no  is  the  ambient  plasma  density.  The 
quiver  current  7,  acts  as  the  source  term  in  the  wave  equation  which  drives  the  harmonic 
radiation,  (V*  —  =  5,,  where 

Sq  kpOq  cos  ^o  (1  +  dg  cos*  0o)  (2) 

and  <ifp  =  ckp  =  (47re*no/m)^/*  is  the  plasma  frequency.  In  the  limit  dj  1,  the  de¬ 
nominator  in  Eq.  (2)  may  be  expanded  and  the  component  driving  the  odd  harmonic 
(a>/v  =  Nuq)  may  be  determined.  For  example,  the  ratio  of  the  third  harmonic  power  to 
the  fundamental  is  P3/P0  {alkpL/16ko)^  =  Rq,  where  L  is  the  laser-plasma  interaction 
length  and  L  <  Ld  has  been  assumed,  where  Ld  is  the  detuning  distance  (discussed  below). 
Hence,  in  the  quiver  limit,  Pj  ~  L^n\Il. 

The  quiver  model  usumed  that  the  plasma  response  is  dominated  by  the  electron 
quiver  motion.  This  is  an  oversimplification  and  collective  plasma  effects,  i.e.,  the  plasma 
density  response,  cannot  be  neglected.  Including  the  density  response,  6n(z,t),  in  the 
transverse  current  gives  7j.  =  -e(no  +  dn)vx.  Letting  t>x  =  v,  gives  S  =  <1*00(1  + dn/no  - 
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®o/2)>  where  the  term  O0/2  arises  from  expanding  the  relativistic  factor,  assuming  aj  1- 
Using  ID  cold  fluid  theory,  the  density  response  Sn  may  be  calculated,  giving 

5  ~  fcjoo  cos^o  ~  ^  ~  cos  2^0^  •  (3) 

Hence,  the  effect  of  including  the  density  response  is  to  reduce  [7],  [8]  the  source  term 
for  the  third  harmonic,  ~  exp(±3i^o)»  hy  the  factor  3a;*/4a>2  -C  1  as  compared  to  the 
quiver  model.  The  power  in  the  third  harmonic  will  be  reduced  by  the  square  of  this 
factor,  P3/P0  =  Rq{2u)^/4tt>Q^)^ .  Hence,  in  the  ID  limit,  P3  ~  L^n^Iq.  Physically,  this 
reduction  arises  from  the  longitudinal  ponderomotive  force,  Fg  ~  —  ~  l;odosin29o- 

This  modulates  the  density,  Snfno  ~  exp(±2t^o)»  in  such  a  way  that  it  nearly  cancels  (to 
•  order  the  contribution  from  the  relativistic  factor,  al/2.  In  the  absence  of  the  space 

charge  potential,  Up  —*  0,  this  cancellation  is  exact  and  no  third  harmonics  are  generated 
in  the  ID  limit. 

The  harmonic  radiation  will  reach  a  maximum  amplitude  after  a  detuning  distance 
[1],  [19],  [20],  L  =  Li.  The  phase  velocity  Vp  of  an  electromagnetic  wave  of  frequency  ui 
in  a  plasma  is  given  by  Vp/c  ~  1  +  u)p/2u)^,  where  ufp/u>^  <  1  and  oj  <C  1  have  been 
assumed.  The  phaise  velocity  of  the  incident  laser,  =  a^o,  is  greater  than  that  of  the  third 
harmonic,  u>  =  3ci;o.  Hence,  the  third  harmonic,  which  is  being  driven  by  the  incident 
laser,  eventually  becomes  out  of  phase  with  the  incident  laser.  The  maximum  amplitude 
of  the  third  harmonic  occurs  after  the  detuning  distance  defined  by  Ld^Vpfc  =  A3/2,  where 
Avj,  =  AcUp/9ul  is  the  difference  in  the  phase  velocities  of  the  incident  and  third  harmonic 
fields  and  As  ~  Ao/3  is  the  wavelength  of  the  third  harmonic.  Hence,  La  ~  3A*/8Ao,  where 
Ap  =  2irfkp  is  the  plasma  wavelength.  The  maximum  amplitude  of  the  third  harmonic 
power  at  saturation  may  be  estimated  by  setting  L  =  L^,  giving  P3/P0  —  (^Ao/Ap)^  %. 
Hence,  P3  ~  nllq.  As  an  example,  consider  a  Aq  =  1  pm  laser  with  /o  =  5  x  10^^  (do  ~  0.6) 
and  a  plasma  of  density  no  =  10^*  cm~’  (Ap  ~  11  pm).  The  third  harmonic  power  is  given 
by  P3/P0  —  9  ^  10“*.  The  detuning  length  is  prohibitively  short,  Li  ~  45  pm.  The  third 
harmonic  pulse  length  is  approximately  equal  to  that  of  the  incident  laser  pulse. 
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II.  NONLINEAR  FORMULATION 


The  ID  i’.elds  associated  with  the  laser-plasma  interaction  can  be  described  by  the 
normalized  transverse  vector  potential,  a(2,t)  =  cAj.(-s,<)/mc*,  and  the  normalized  scaler 
potential,  ^(z,t)  =  e#(2,t)/mc^.  Coulomb  gauge  will  be  used,  V  •  A  =  0,  which  implies 
A,  =  0.  It  is  convenient  to  introduce  the  independent  variables  ^  =  z  —  cfift  and  t  =  t, 
where  fit  is  the  normalized  transform  velocity  which  will  be  specified  below.  Using  the  t 
variables,  the  normalized  potentials  a  and  ^  satisfy  the  ID  wave  equation  and  Poisson’s 
equation,  which  are  given  by 


7i  c  didr  dr^ )  ^  ^7no*’ 


where  7^“*  =  1  —  fif ,  Up  =  ckp  =  (4irc^no/Trt)^'^*  is  the  ambient  plasma  frequency,  n(^,r) 
is  the  plasma  electron  density,  no  is  the  ambient  density,  7(^,t)  =  (1  —  is  the 

relativistic  factor  associated  with  the  plasma  electrons  and  ;9(^,t)  =  v/c  is  the  normalized 
electron  fluid  velocity.  In  deriving  the  right  side  of  Eq.  (4),  use  was  made  of  the  fact  that 
the  transverse  canonical  momentum  is  invariant  and  that  prior  to  the  laser  interaction 
the  plasma  is  assumed  to  be  stationary,  i.e.,  fix  —  o/t-  Also,  the  ions  are  assumed  to  be 
stationary. 

The  electron  fluid  quantities  n,  fiz  and  7  are  assumed  to  satisfy  the  cold,  relativistic 
fluid  equations.  Using  the  ^,r  variables,  the  continuity  equation  and  the  axial  momentum 
equation  may  be  written  as 

=  (6) 

^|7(l-AA)-«>]  =  -i|;(7|9.),  (7) 

respectively.  Equations  (4)-(7),  together  with  7  =  (H-o*)*/*/(l  —fi\Y^^,  form  a  complete 
set  of  fully  nonlinear,  relativistic  cold  fluid  equations  which  describe  the  ID  laser-plasma 
interaction.  The  ID  assumption  is  valid  provided  that  the  radiation  spot  size  t,  is  large 
compared  to  the  plasma  wavelength  Ap  =  2w/kp,  i.e.,  r,  Ap.  The  cold  fluid  assumption 
is  valid  provided  that  (i)  the  electron  quiver  velocity  is  much  greater  than  the  transverse 
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thermal  velocity  and  (ii)  the  phase  velocity  of  the  driven  plasma  oscillations  is  much  greater 
than  the  axial  thermal  velocity,  as  is  the  case  in  the  following. 

The  primary  focus  of  this  Section  is  to  determine  the  coherent  harmonic  radiation  ex¬ 
cited  by  a  long  pulse  incident  laser  field,  cti  »  Xp.  To  study  coherent  harmonic  radiation, 
it  is  convenient  to  set  the  transform  velocity,  /dt,  equal  to  the  phase  velocity,  of  the 
incident  laser  field,  fit  =  fipk-  incident  laser  field  is  of  the  form  oq  =  do  cos  where 
kfii  =  ko{z  —  cfipitt),  fiph  =  «o/cfco  is  the  normalized  phase  velocity,  wo  is  the  frequency  and 
A;o  is  the  wavenumber  of  the  pump  laser  field.  The  pump  laser  amplitude  is  slowly  varying, 
|ddo/5^|  ~  is  assumed  to  be  independent  of  t,  i.e.,  pump  depletion  effects  are 

neglected.  The  effects  of  diffraction  are  consider  in  Section  III.  In  determining  the  plasma 
response  to  an  incident  laser  field  of  the  form  ao  =  ao(0»  *be  r  derivatives  are  dropped 
in  Eqs.  (6)-(7).  Strictly  speaking,  for  a  general  laser  field  of  the  form  oo  =  ao{^,T),  the  t 
derivatives  may  be  neglected  provided  that  the  transit  time  of  the  electrons  through  the 
laser  pulse,  (=  laser  pulse  duration),  is  small  compared  to  the  evolution  time,  te,  of  the 
laser  envelope,  ti  C  te,  i.e.,  the  quasi-static  approximation  [8].  Assuming  a  quasi-static 
plasma  response,  Eqs.  (6)-(7)  imply  the  existence  of  two  constants  of  the  motion. 


n{fit  -  fix)  =  fitno,  (8a) 

^{i-fitfix)-<i>  =  i^  m 

Equations  (8a)-(8b)  may  be  solved  to  give  expressions  for  the  fluid  quantities  n,  fi^  and  7 
in  terms  of  the  fields  a  and  4>, 


n/no  =  nflfit  [(1  -  -  ^t]  , 

7  =  7?(1  +  ^)  [l-W-l/V’S?)^/*], 
n/7no  =  fit{l+4>)~^il  - 


(9a) 

(96) 

(9c) 

(9d) 


where  =  (1  -f  ^)*/(l  +  o’).  Notice  that  for  the  case  fit  =  fiph,  7«  ^  =  Iph  < 
0ih>^  for  a  radiation  field  in  a  plasma. 
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Using  the  above  expressions  for  the  fluid  quantities,  expressions  for  the  normalized 
transverse  plasma  current,  fcpna/no7,  and  the  normalized  charge  perturbation,  kp{n/no  — 
1),  can  be  found  emd  inserted  into  the  wave  equation  smd  Poisson’s  equation,  Eqs.  (4)-(5). 
This  results  in  two  coupled  nonlinear  equations  for  the  potentials  a  and  (f>, 

1 


c  didr  c^dr^ 


[/J.(i  -  i/V-S?)-'’  - 1]  • 


1 

(10) 

-']• 

(11) 

Equations  (lO)'(ll)  completely  describe  the  nonlinear  generation  of  relativistic  harmonic 
radiation  in  ID  within  the  quasi-static  approximation,  i.e.,  te-  For  the  case  where 

the  transform  velocity,  is  set  equal  to  the  phase  velocity,  of  the  incident  laser,  Eqs. 
(lO)-(ll)  may  be  used  to  analyze  relativistic  harmonic  generation,  as  is  done  in  detail  in 
the  following.  For  the  case  where  0t  is  set  equal  to  the  group  velocity,  ,  of  a  laser  pulse 
with  a  pulse  length  ^  \p,  Eq.  (11)  may  be  used  to  analyze  laser  wakefield  generation 
[5]- [9].  The  limiting  case  o{  dt  =  1  has  been  analyzed  in  detail  in  Ref.  [8]. 

It  is  of  interest  to  define  the  nonlinear  index  of  refraction  ur  by  setting  the  right  side 
of  the  wave  equation,  Eq.  (10),  equal  to  (a>’/c*)(l  —  n^)a,  which  gives 


"*  “  a.»(i + ^)  V  thf  se )  ■ 


(12) 


In  particular,  the  slow  part  of  ur  determines  the  dispersion  relation  for  the  radiation 
field  a,  whereas  the  fast  part  of  ur  determines  the  generation  of  harmonic  radiation,  sis  is 
discussed  below. 


A.  Dispersion  Relation 

To  describe  the  generation  of  relativistic  harmonics,  the  electrostatic  potential  <f>  is 
separated  into  slow  and  fast  components,  <!>  =  +^/,  where  (f>,  is  approximately  constant, 

\d(l>tldi\  and  \d<f>fldi\  ~  The  fast  component  <f>f  contains  the  harmonic 

content  and  is  important  in  determining  the  source  current  which  drives  the  harmonic 
radiation.  The  slow  component  <f>,  is  important  in  determining  the  dispersive  character¬ 
istics  of  the  harmonic  radiation  as  well  as  of  the  pump  laser  field.  In  the  following,  it  is 
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assumed  that  ^*/(l  +  <C  1.  Furthermore,  it  is  assumed  that  the  vector  potential  of 
the  harmonic  radiation  a^v,  where  N  is  the  harmonic  number,  is  much  less  than  that  of 
the  pump  laser,  i,e.,  |ojv|/|do|  <  1. 

The  dispersion  relation  is  determined  by  examining  the  slow  part  of  the  Eq.  (12). 
Assuming  ^  1>  the  dispersion  relation  for  the  radiation  field,  a,  to 

leading  order  (i.e.,  neglecting  is  given  by 

=  +  (13) 

where  a  =  aexp(tX;^2  —  has  oeen  assumed  with  a  slowly  varying  compared  to 

the  phase  factor  cxp{ikf,z  —  Equation  (13)  holds  for  the  pump  laser  field  (u;,,  = 

=  ^o)  M  well  as  for  the  various  harmonics.  Notice  that  the  transform  velocity  is 
the  phase  velocity  of  the  incident  laser  field,  fipk  =  u>o/cko,  and,  hence,  7“^  =  1  -  (3^^  = 
— fcp/Ao(l  +  ^«)  <  0-  Furthermore,  for  a  nonevolving,  long  incident  laser  pulse  of  the  form 
Uo  =  do  exp(tX;o'2!  —  ‘‘^ot),  the  quantity  (1  +  </>,)  is  approximately  constant,  as  is  discussed 
in  the  following. 

Consider  a  long  pump  laser  pulse  with  a  slowly  varying  envelope,  |5do/d^|  ~  (do/cr£,|, 
where  is  the  laser  pulse  duration,  propagating  in  a  plasma  with  cr^,  >  Ap,Ao.  The 
slow  part  of  the  potential  may  be  determined  from  examining  the  slow  part  of  Poisson’s 
equation,  Eq.  (11)>  In  the  long  pulse  limit,  d^<f>t/d(^  may  be  neglected  and  Eq.  (11)  implies 
that,  to  leading  order  (i.e.,  neglecting  ^/)  [8],  [21], 

(1  +  W  =  (1+W'’-  (H) 

In  particular,  notice  that  in  the  long  pulse  limit,  cti  Ap,Ao,  propagation  (i.e.,  real  k) 
requires  u»p/a>*(l  +  <f>i)  <  1.  This  implies  that  for  intense  pump  laser  fields  with  do  >  1, 
propagation  in  an  overdense  plasma  in  which  u>p/u;^  >  1  may  be  possible.  Physically, 
7,  =  1  +  =  (1  +  al/2y^^,  and  the  reduction  in  the  effective  plasma  frequency, 

is  due  to  the  relativistic  quiver  motion  of  the  electrons. 

It  should  be  emphasized  that  Eq.  (14)  only  applies  to  the  long  pulse  limit  ct^  > 
Ap,  Ao.  For  Ti  =  1  psec,  ti  »  Ap  implies  plasma  densities  no  '>  10**  cm“®.  In  the 
short  pulse  limit,  Ap  >  ctl  >  Ao,  it  can  be  shown  [8],  [21]  that  |^,|  -C  1  provided 
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ctl  C  Ap/(1  +  0^/2)^/^.  Physically,  \4>t\  <  1  holds  for  pulses  sufficiently  shorter  than  \p, 
since  \p  is  the  characteristic  length  scale  (in  the  ^  frame)  for  collective  electron  motion 
(i.e.,  collective  electron  motion  leading  to  charge  separation  and  significant  values  of  <f>, 
does  not  occur  on  time  scales  sufficiently  shorter  than  l/up). 

B.  Plasma  Response 

The  source  current  driving  the  harmonic  radiation  is  determined  by  the  harmonic 
content  of  tir.  As  indicated  by  Eq.  (12),  riR  may  be  specified  entirely  in  terms  of  (j). 
Hence,  it  is  necessary  to  determine  the  various  harmonic  components  of  <t>.  This  is  done 
by  analyzing  the  plasma  response  to  the  pump  laser  field,  oq  =  do  cos  via  Poisson’s 
equation,  Eq.  (11).  The  various  harmonic  components  of  where  4>2i  ~ 

cos(2^fcoO  ^  =  1,2,3...,  may  be  determined  analytically  in  the  limit  c  =  kp/kl{l  + 
<t>g)  <  1.  In  this  limit  \<f>2tl{l  +  ^j)l  ~  c^,  and  the  various  harmonic  components  o{  <f>f  may 
be  solved  order  by  order  in  e. 

To  solve  for  via  Eq.  (11),  two  expansions  will  be  used.  Assuming  |(V’7pfc)~*l  ^ 


gives 


^(2£-l)l!<^  1  \ 


(15a) 


where  (2£  -  1)!!  =  1  •  3  • ...  •  (2£  -  1)  and  (2£)!!  =  2  •  4  • •  (2£).  Assuming  |^//(1  +  <f>,)\  <  e. 


gives 


(1  + 
7jh(l  + 


^  2£(2£+l)...(2£  +  Tn-  1)  -<f>f 

^  m!  .(l  +  <^.) 


(156) 


Letting  0/  =  where  <l>2t  ~  c^co8(2£fcoO  for  ^  =  1>2,3...,  and  using  the  expansions 

in  Eqs.  (15a)  and  (15b),  allows  <f>2t  to  be  solved  order  by  order  in  c.  For  example,  <f>2  (order 
e),  <t>4  (order  c*)  and  (order  e’)  are  given  by 


_  kl{al)2 

2{\  +  4>.y' 


(16a) 


^2(00)2 


M  (l  +  <?^.)’  87j,(l+<^.)^J  ’ 


(166) 
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L  _  3.^1  ,  3(ag)2  /  _  5{alh  \ 

(1  P  2il  +  <f>,)  ^  2-rl^{l  +  <i>,y  247*,(1 


(16c) 


where  the  subscript  (Q)<  refers  to  the  harmonic  component  of  Q,  i.e.,  {Q)t  ~  exp(ffcoO- 
These  equations  can  be  solved  iteratively  to  yield 


kpCil  cos2fco^ 

l6fc2(TT^’ 


<f>i  = 

<t>6  = 


fepdfl  cos4A:o^ 

2»fc*(i + <i>,r  ’ 

ITfcpd®  cos  6A:o^ 


(17a) 

(176) 

(17c) 


For  order  of  magnitude  scaling  purposes,  Eqs.  (17a)-(17c)  imply  that  the  2^*^  harmonic 
component  ol  4>f  scales  as 


,  ,  , ,, (2£  -  1)!!  f  al^  cos 2ikoi 

1)  ^2£)!!  V2jbo/  2£2(l  +  ^,)3^-i- 


(18) 


Knowing  the  harmonic  components  of  it  is  possible  to  determine  the  harmonic 
components  of  the  source  current  5  which  drives  the  harmonic  radiation.  The  source 
current  S  is  given  by  the  right  side  of  Eq.  (10), 


5  = 


1 


(19) 


Using  the  expansion 


(1+^)-'  =  (1+#.)-'  Y.  (2») 

m=0 


allows  the  harmonic  components  of  S  to  be  solved  order  by  order  in  e.  For  example,  the 
third,  5s,  fifth,  5$,  and  seventh,  57,  harmonic  components  are  given  by 


53  = 


fepg 


<i>2  ^  \  d^<i>2 


(1+^.)  [  (1  +  ^.)  klyl^  dC  J 


5s  = 


-kla 

^4 

<f>2  ^ 

(  ^^4 

h  d^<l>2\ 

(1  +  ^») 

(1  +  <f>t) 

(1 + <i>,y 

{ee 

(1  +  ^.)  de ) 

(21a) 

(216) 
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57  = 


1 


+ 


<i>l 


2(^204 


d^<f>i  <f>2  d^<p4  f  4>4 


khlk[d(^  {i  +  <i>.)de  \{1 +  <(>.) 


P  'ph, 

Using  Eqs.  (17a)-(17c)  give 


-G 


<i>\ 

0-  +  <t>.? 


Si 

5s 


3k*al  cosSko^ 
“32fc2(l+<^,)4’ 
ZkpCil  cosSfeo^ 


5k*al cos  7ko^ 

3‘2^^kl{l  +  <f>,y^- 


(21c) 


(22a) 

(226) 

(22c) 


For  order  of  magnitude  scaling  purposes,  Eqs.  (22a)-(22c)  imply  that  the  {2£  +  1)‘^  har¬ 
monic  component  of  5  scales  as 


S2i+i  ~  (-1)'*; 


(2£-l)!! 

(2£)!! 


/  kp\^^  cos(2£  -I-  l)fco^ 
\2ko)  + 


(23) 


C.  Harmonic  Radiation 


The  above  expressions  for  the  harmonic  components  of  the  source  current,  S2/+1,  may 
be  used  in  the  wave  equation,  Eq.  (10),  to  determine  the  growth  and  saturation  of  the 
relativistic  harmonic  radiation.  It  is  convenient  to  represent  the  harmonic  radiation 
field,  aj^ ,  by  the  form 

(24) 


with  Aufjf  =  Uff  —  ^pf^ckff  and  k^(  —  =  k^z  —  where  k^  and  u,,  are  the 

wavenumber  and  frequency  of  the  N*^  harmonic  radiation  field.  Inserting  this  into  Eq. 
(10)  gives  a  reduced  wave  equation  for  the  harmonic  amplitude  a,^. 


d  . 


=  -c^5,,exp(tAu;,^T), 


(25) 


where  (iff,  and  k„  satisfy  the  dispersion  relation,  Eq.  (13),  with  k,,  =  Nko  and  N  = 
2^  -f-  1.  The  amplitude  of  the  source  current  for  the  N*^  harmonic,  Sf,,  is  given  by  S,,  = 
Sf,  exp(tArjboO»  S„  is  given  in  Eqs.  (22)-(23). 
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The  reduced  wave  equation,  Eq.  (25),  has  the  solution 


2c  [e3cp(iAa;^T)  -  1] 


where 


a - - 


For  small  times  |Awj^t|  •<  1,  the  harmonics  grow  linearly  in  time  [8], 

The  N*'^  harmonic  reaches  its  maximum  amplitude  [1],  [19],  [20]  after  a  dephasing 
length  cra^  =  7rc/|Act»j,  |,  where  Aofj^  =  -  ^ph)  and  =  (v^/ck^,  i.e,, 

V  k^  /  k^ 

Au;^  =  Ncko  1  +  ■■■  ,  -  1  +  .  (29) 

\  Nniii+4>s)  V  kl{i  +  4>.) 


For  A;^/l!o(l  +  ^*)  ■C  1,  Eq.  (29)  may  be  expanded  giving 

The  maximum  amplitude  of  the  harmonic  at  cr  ~  crj^  is  given  by 

_  215,1(1  + d.) 

Klm«.  -  *,2(JV2_1)  •  V31) 

In  particular,  for  the  third,  fifth  and  seventh  harmonics, 

1®’'"*“*“  2^jb2(l  +  <^.)3,  ^  ’ 

"  9-2i3jb»(l+0,)'»’ 

An  order  of  magnitude  expression  for  the  maximum  amplitude  of  the  (2£  +1)  harmonic  is 


(326) 


(32c) 


given  by 


2{2t-mkj2hY^af^^ 

|a2/+jU«*  -  ^(2£)!f  [(2£  +  1)2  -  l](l  +  <^2)3/- 
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Equations  (32)-(33)  are  valid  for  all  do,  including  do  >  1.  It  is  easily  verified  that 
l®jvlmox  ’C  do.  A  plot  of  the  harmonic  radiation  amplitude  (d^^  |  verses  the  normalized 
propagation  distance  ct/Aq  is  shown  in  Fig.  2  for  the  first  three  harmonics  for  the  pa¬ 
rameters  Ao  =  1  fixa,  Ap  =  3.4  fim  (no  =  10*”  cm"*)  and  do  =  2.2  (Jo  =  6.7  x  10** 
W/cm*). 

At  saturation  for  the  harmonic,  i.e.,  cr  =  ct^^,  the  ratio  of  the  power  in  the 
N*^  harmonic,  P^,  to  the  power  in  the  pump  radiation  field,  Po,  is  given  by  P^/Po  = 
Using  Eqs.  (32)-(33),  this  can  be  written  as 


Pn 

Po 


2(N-1) 


H) 


(34) 


where  Cn  m*  constants  which  decrease  rapidly  with  increasing  harmonic  number,  i.e., 
C7s  =  4.9  X  10  Cs  =  2.4  x  10“®  and  C^  =  2.3  x  10"^.  Furthermore,  driven  relativistic 
harmonic  generation  is  a  nonresonant  interaction;  hence,  the  process  is  not  sensitive  to 
thermal  plasma  effects. 

As  an  example,  consider  an  incident  laser  with  Po  =  10  TW,  Ao  =  1  pm  and  spot  size 
ro  =  10  pm  (do  =  2.2),  interacting  with  a  plasma  of  density  no  =  10*”  cm~*  (Ap  =  3.4 
pm).  For  the  third  (and  fifth)  harmonic,  P/v/Po  =  2.2  x  10~®  (4.6  x  10~*”)  and  La  =  8.0 
pm  (4.3  pm).  Hence,  220  MW  (4.6  kW)  should  be  observed  at  a  wavelength  of  A  =  3300 
A  (2000  A).  Clearly,  the  limitations  due  to  phase  detuning  are  restrictive.  If  a  scheme  for 
phase  matching  could  be  conceived,  the  interaction  distance,  L,  and  thus  the  harmonic 
power,  P/v  ~  L^,  could  be  increased. 
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III.  DIFFRACTIVE  EFFECTS 


The  above  nonlinear  theory  indicates  that  the  N*''  harmonic  saturates  (i.e.,  reaches 
maximum  amplitude)  after  a  detuning  distance  cr^f,  =  Trc/jAwj^l,  where  =  ck^{f3f,  — 
fiph)-  In  particular,  for  the  third  harmonic,  ~  3Ap(l  +  ^,)/8Ao.  Physically,  saturation 
by  detuning  arises  due  to  the  fact  that  the  phase  velocity  of  the  fundamental  pump  radi¬ 
ation  is  greater  that  than  the  phase  velocity  of  the  harmonic  radiation,  i.e.,  I3p\  >  0^-  In 
a  realistic  3D  configuration,  however,  diffractive  effects  will  limit  the  effective  radiation- 
plasma  interaction  length  to  a  few  vacuum  Rayleigh  lengths,  Zr,  where  Zr  =  korl/2 
and  To  is  the  minimum  spot  size  of  the  radiation  field  (assumed  to  be  Gaussian).  For 
interaction  distances  L  i^Zr,  diffraction  effects  are  unimportant  and  the  ID  theory  is 
an  adequate  description.  In  the  limit  itZr  <  L,  however,  and  in  particular  itZr  <  cTtif, , 
diffractive  effects  are  important  and  must  be  included  in  the  analysis. 

In  the  following,  the  generation  of  second  and  third  harmonic  radiation  is  analyzed 
using  3D  relativistic  fluid  equations  in  the  limit  Oq  ‘C  1.  Higher  order  harmonic  generation 
can  be  consider  by  solving  the  fluid  equations  to  higher  order  in  a*.  Furthermore,  the 
radiation  fields  are  assumed  to  undergo  vacuum  diffraction.  Notice  that  in  the  regime 
in  which  3D  effects  are  important,  the  condition  itZr  <  CTd^  implies  S  =  2k^rl/Z  <  1. 
Hence,  the  effects  of  relativistic  optical  guiding  [3]-[10],  which  become  important  when 
P/Pc  =  k^r\a\/Z2'^\,  may  be  neglected  in  the  3D  regime  when  1.  Assuming 

vacuum  diffraction,  a  long  pulse,  Gaussian,  incident  radiation  field  evolves  according  to 
o(r,z,<)  =  o,  expt^o*  where  =  k^z  —  u>ot  and  the  radiation  envelope  is  given  by  [24] 

o,(r,z)  =  ^^exp  -(1  -  to)^  -  itan"^  a  ,  (35) 

where  r,  =  ro(l  +  is  the  radiation  spot  size,  a  =  z/Zr  and  z  =  0  is  the  location  of 

the  laser  focus  at  which  r,  =  tq. 

A.  3D  Formulation 

The  3D  pump  laser  interaction  with  the  plasma  electrons  will  be  modelled  using  the 
cold,  relativistic  fluid  equations.  In  particular,  the  momentum  equation  and  the  continuity 
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(36) 


equation  may  be  written  as 

-_u  =  V^+--.--ux(Vxa), 

+  ^  •(?“)  =  (•.  (37) 

respectively,  where  u  =  p/mc  is  the  normalized  momentum,  7  =  (1  is  the  relativis¬ 

tic  factor,  and  p  =  n/(no7).  The  evolution  of  the  normalized  potentials  are  determined 
from  the  wave  equation  and  Poisson’s  equation,  which  may  be  written  as 


W  =  -  1), 


respectively,  where  Coulomb  gauge,  V  •  a  =  0,  has  been  assumed. 

To  study  the  generation  of  the  third  harmonic,  the  above  equations  will  be  solved 
perturbatively,  i.e.,  order  by  order  in  the  pump  laser  amplitude  do,  assuming  dg  <C  1. 
The  various  plasma  quantities,  denoted  by  Q,  will  be  represented  by  an  expansion  Q  = 
Qt>  +  Q\  +  Q2  where  Qn  ~  aj .  The  pump  laser  vector  potential  is  assumed  to  be 
known  and  given  by  Eojl  =  o,  expt<>o«*>  where  a„(r,  z)  is  the  vacuum  solution  given  by 
Eq.  (35).  The  effects  of  pump  depletion  and  of  various  laser-plasma  instabilities  will  be 
neglected.  The  zeroth  order  plasma  quantities  (in  the  absence  of  the  pump  laser)  are  given 
by  Uo  =  0,  7o  =  1,  Po  =  1  *nd  ^0  =  the  plasma  is  assumed  to  be  initially  uniform, 

stationary  and  neutral.  To  first  order  in  the  normalized  vector  potential,  Eq.  (36)  implies 
that  Ui  =  Ho,  which  is  simply  the  quiver  motion  of  the  electrons.  Furthermore,  Eqs.  (36)- 
(39)  imply  Pi  =  7i  =  =  0.  The  first  order  form  of  the  wave  equation,  Eq.  (38),  implies 

that  the  pump  laser  field  obeys  the  dispersion  relation  Wg/c*  =  fcg  -{-  fcp. 

To  second  order,  Eqs.  (36)*(39)  imply  that  quantities  ^3,  U2  and  pz  are  related  to  Og 


by 


(40a) 

(406) 

(40c) 
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Furthermore,  Eqs.  (40a)  and  (40c)  imply  that  the  second  order  wave  equation  describing 
the  generation  of  the  second  harmonic  may  be  written  as 


c2  dt^ 


-  hpj  »2  =  0. 


(41) 


Hence,  no  second  harmonic  radiation  will  be  generated. 

To  third  order,  Eqs.  (36)-(39)  imply  that  quantities  ^3  and  U3  are  related  to  a®  by 


vV,  =  * 


V^(ao  •  U2)  -  •  (^0^2) 


1  ^  r,/ 

-^Us  =  V^3- V(ao-U2), 


(42a) 

(42ft) 


Furthermore,  the  third  order  transverse  wave  equation,  describing  the  generation  of  third 
harmonic  radiation,  is  given  by 


~  ^  (*Jp2ao±  +  fcju3±  +  ,  (43) 

where  (...)n  signifies  the  n*^  harmonic  component. 

In  particular,  to  determine  the  source  term  driving  the  third  harmonic  radiation,  it  is 
necessary  to  determine  the  fast  (i.e.,  third  harmonic)  components  of  ^3  and  Us.  Assuming 
^3iU3  exp3t0O)  Eqs.  (42a)  and  (42b)  imply 

u;2 

</>3  ^  •  ^22)3*  (44a) 

ic 

U3  ~  -- — V(ao  •  U2)3,  (44ft) 

where  terms  of  order  k^/kg  •<  1  and  \lr\kl  C  1  have  been  neglected.  It  can  be  shown 
that  the  leading  order  contributions  of  each  of  the  three  terms  on  the  left  side  of  Eq.  (43) 
are  of  the  same  order.  However,  to  leading  order,  k^Usj,  +  Vj.5^3/5ct  =  0,  and,  hence, 
the  wave  equation  describing  third  harmonic  generation  is  given  by 

(v’  -  ~  ajj,  =  tJ(^aox)2,  (45) 

where  higher  order  terms  (of  order  kp/kg  and  1/r^k^)  have  been  neglected.  The  transverse 
wave  equation,  Eq.  (45),  along  with  the  second  harmonic  component  of  pj, 
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completely  determine  the  generation  of  third  harmonic  radiation  in  3D  in  the  limit  <C  1 . 

B.  Third  Harmonic  Generation 


The  pump  laser  is  assumed  to  be  a  long  pulse,  Gaussian  laser  beam  which  is  diffracting 
according  to  the  relation  ao(r,  z,  t)  =  o,  expi^Ot  where  o„(r,  z)  is  the  envelope  given  by  Eq. 
(35)  and  u;^/c’  =  4-  As*.  Using  this  form  in  Eq.  (46)  gives 


2rlkl 


(47) 


Hence,  the  source  term  driving  the  third  harmonic  radiation  is  given  by 


*:J(p2<»o)s  - 


ife* 


3  2t2  ,  1 

+ 


(1-f  ta)J 


al  exp  St^o  • 


(48) 


It  is  convenient  to  denote  the  third  harmonic  radiation  field,  avs{r,z,t),  by  the  fol¬ 
lowing  form. 


®8(r,z,<)  =  /(2)o,8(r,z)expids, 


0,3(2,  r)  =  —  exp 


— (1  —  »as)-5 - »tan“^  ots 


•s 


(49a) 

(496) 


where  0$  =  k^z  —  ufsi,  u>s  and  ks  are  the  third  harmonic  frequency  and  wavenumber, 
r»8  =  ra(l  +  is  the  third  harmonic  spot  size,  03  =  2/^113,  Zrs  =  k3r$/2  is  the 

Rayleigh  length  associated  with  the  third  harmonic,  and  u\/c^  ■=  k\  ■\-k'^.  Inserting  this 
into  the  wave  equation,  Elq.  (45),  gives 


Q 

'2tfcso,3  exp(itfs)^/  *J(p20o)3,  (50) 

where  the  right  side  is  given  by  Eq.  (48).  Letting  u;3  =  3<4;o  and  r|  =r\/Z  implies  Zm  ~  Zr 
and 

Bz’  -  24rjt’ 
where  Afc  =  3fco  —  fcs  =  — 4fcJ/3fco. 

For  a  pump  laser  pulse  which  interacts  with  a  uniform  plasma  extending  from  Zmin  < 
*  <  *ma*»  Eq.  (51)  can  be  integrated  giving 


exp(tAA;2) 
(1  +  m)»  ’ 


(51) 


)  -  Otmin)]  9 


(52o) 


h{8, 


(526) 


(1  +  ta)  y  (1  +  ta) 

where  \f\D\  —  \Za\k^/2^ kl\  is  the  maximum  amplitude  of  the  third  harmonic  in  the  ID 
limit,  6  =  \^k\Zll  =  2kpVQ/3  (6  =  —  ^min/ O^max  —  ^max/^R  and  Ei 

is  the  exponential  integrsd. 

In  particular,  consider  a  finite  slab  plasma  centered  about  the  focal  point  of  the  pump 
laser  (chosen  to  be  x  =  0)  which  extends  from  z^in  =  -*o  to  z^ax  =  2o-  The  amplitude 
function  /  of  the  third  harmonic  signal  emerging  from  the  slab  plasma  at  z  =  zo  Is  given 
by 

/  =  2|/iD|Im[/i(6,ao)],  (53) 


where  Im[h(6,ao)]  is  the  imaginary  part  of  h{S,ao)  and  oq  =  zqIZr.  The  ratio  of  the 
third  harmonic  power  to  the  pump  laser  power  is  given  by 


P,IP«  =  =  3|/lViS-  (64) 

Hence,  for  a  slab  plasma  centered  about  the  laser  focus,  Ps  ~  4|Im  [h(6,  ao)]  1*.  A  plot  of 
the  function  J5r(^,ao)  =  4|Im[h(6,ao)]  ^  verses  the  normalized  plasma  slab  width,  qo,  is 
shown  in  Fig.  3  for  5  =  0.1,  1.0  and  10.0.  Notice  that  for  6  =  10.0  (cr^,  «C  X’Zfl)  and 
ao  C  1,  the  first  peak,  which  occurs  when  the  slab  width  is  equal  to  the  detuning  length, 
2zo  =  cTiy,  is  close  to  the  ID  value  of  H{6,ao)  =  1.  For  8  =  0.1  (era,  >  ’tZk),  the 
maximum  amplitude  of  the  third  harmonic  occurs  for  a  plasma  length  of  2:^  Zjt,  at  which 
H^(^,ao)  ^  1/3.  In  the  limit  of  a  slab  plasma  with  a  spatial  extent  large  compared  to 
Zr,  i.e.,  zo  >  Zr  (ao  »!),/—»  //,  where  //  =  |/ix?|a^*  sin(6ao).  Hence,  //  — ►  0  as 
ao  — »  00,  i.e.,  the  third  harmonic  radiation  emerging  from  a  plasma  centered  about  the 
laser  focus  vanishes  for  plasma  dimensions  large  compared  to  Zr.  The  conclusion  that 
no  third  harmonic  is  generated  by  an  infinite  medium  is  also  the  case  for  third  harmonic 
generation  due  to  the  nonlinear  susceptibility  associated  with  bound  electrons  [25]. 

A  finite  third  harmonic  signal,  however,  may  be  detected  by  focusing  the  pump  laser 
on  the  trailing  edge  (z  =  0)  of  a  plasma  slab  extending  from  Zmin  <  ^  ^  Zmax  =  0,  with 
dimensions  large  compared  to  Zr,  i.e.,  z„,j„  C  —Zr.  In  practice,  this  may  be  achieved  by 
focusing  the  pump  laser  pulse  on  the  trailing  edge  of  a  pre-ionized  gas  get.  Approximating 
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(55a) 


Omin  =  -oo  and  a^o*  =  0  gives  f  =  fs,  where 

giS)  =  (S/8)[l-Se‘E,(S)].  (555) 

Furthermore,  it  can  be  shown  that  0  <  1  —  Se^Ei{S)  <  (1  +  5)~^'  Hence,  -C  1  and  may 
be  neglected.  The  amplitude  of  the  third  harmonic  emerging  from  the  plasma  is  given  by 
|as|  =  |/s||avs|.  At  the  laser  focus,  this  is  a  factor  of  4/9  smaller  than  that  obtained  from 
the  ID  theory  in  the  limit  o*  C  1.  Hence,  for  a  semi-infinite  plasma  including  the  effects 
of  diffraction  in  the  limit  dg  1,  the  ratio  of  third  harmonic  power  to  the  pump  laser 
power,  Pi/Po  =  3|/||/dg,  is  given  by 

Ps/Po  ^  3  X  10-‘‘(Aoao/Ap)^  (56) 

which  is  a  factor  of  (4/9)* /3  2;  1/15  smaller  than  the  corresponding  ID  expression  for  the 
maximum  power. 

C.  Second  Harmonic  Generation 

The  above  results  indicate  that  for  an  initially  uniform  plasma  density,  no  second 
harmonics  are  excited.  This  is  true  even  when  the  effects  of  a  finite  incident  laser  spot  size, 
ro,  and  the  transverse  gradients  associated  with  self-consistent  3D  density  perturbation 
are  included,  i.e.,  Eq.  (41).  Generation  of  even  harmonics  requires  that  the  plasma  have 
initial  transverse  density  gradients  prior  to  the  arrival  of  the  incident  laser,  i.e.,  ^  0, 

where  no(r)  is  the  initial  plasma  density  profile.  Electrons  undergoing  quiver  motion 
in  the  presence  of  a  density  gradient  produce  density  oscillations.  For  a  pump  laser  of 
the  form  ao  =  oo  ^he  continuity  equation,  to  first  order  in  do  'C  1,  implies  a 

density  oscillation  6nx  of  the  form  6ni  ~  (do  sin®o)V*no/fco-  This  produces  a  source 
current  —  ibJ(5ni/no)do  cos^o  ~  8in25o  which  drives  the  second  harmonic  radiation. 
The  amplitude  of  the  second  harmonic  radiation  may  be  estimated  by  approximating  the 
density  gradient  by  V.no  =;  -no/vp,  assuming  rj  »  rj  and  neglecting  the  effects  of 
diffraction.  At  saturation,  i.e.,  after  a  detuning  length  =  2A*/3Ao,  the  amplitude  of 
the  second  harmonic  is  given  by  |d2|  al/SkoVp.  The  ratio  of  second  harmonic  power 
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to  pump  laser  power  is  P2/P0  —  (aoAo/rp)^  %  and,  hence,  P2  ~  /q-  Not  only  do  initial 
density  gradients  lead  to  the  generation  of  even  harmonics  for  linearly  polarized  incident 
lasers,  it  also  implies  that  circularly  polarized  incident  lasers  will  generate  harmonics,  both 
even  and  odd. 

Transverse  plasma  density  gradients  can  occur  when  a  Gaussian  laser  pulse  produces 
intensity-dependent  ionization  of  a  neutral  gas.  For  laser  pulses  with  peak  intensities 
h  >  Is  (where  Is  is  the  saturation  intensity  for  which  the  gas  is  fully  ionized)  ionizing  a 
uniform  gas,  density  gradients  will  exist  in  the  “halo”  region  about  the  laser  focus.  This 
halo  region  is  the  portion  of  the  interaction  region  in  which  the  gas  is  not  fully  ionized  and 
corresponds  approximately  to  focal  regions  in  which  the  laser  intensity  lies  within  the  band 
Imin  <  I  <  Is  (where  Imin  is  the  minimum  intensity  required  to  produce  ionization).  As 
the  peak  intensity,  /q,  increases,  one  can  show  that  the  volume  of  this  halo  region  increases 
as  /q  ,  where  a  Gaussism  laser  pulse  of  the  form  given  by  Eq.  (35)  has  been  assumed 
(i.e.,  a  spherical  lens  and  a  double  cone  focal  geometry)  [26].  Second  harmonics  will  be 
produced  from  the  halo  region.  The  intensity  scaling  of  the  second  harmonic  power  with 
peak  intensity  will  tend  to  be  dominated  by  effects  associated  with  the  increasing  volume 
of  the  halo  region.  Hence,  for  /q  >  /s>  fh*  volume  effect  implies  a  second  harmonic  power 
scaling  of  P2  /g'  .  Furthermore,  in  the  partially  ionized  halo  region,  the  production  of 
harmonic  radiation  can  be  significantly  enhanced  and/or  dominated  by  atomic  [13]  and 
ionization  [14]  processes  [17],  [18].  In  long-pulse  laaer-plaama  experiments,  filamentation 
may  be  a  dominant  mechanism  in  second  harmonic  generation  [15],  [27]. 


20 


IV.  CONCLUSION 


A  nonlinear  cold  fluid  model  valid  in  the  regime  >  \,  given  by  Eqs.  (10)  and 
(11),  has  been  formulated  and  used  to  analyze  relativistic  harmonic  generation.  The  self- 
consistent  collective  plasma  response  is  included  and  shown  to  significantly  reduce  the 
source  current  driving  the  harmonic  radiation.  A  general  expression  for  the  nonlinear 
index  of  refraction  was  derived,  Eq.  (12),  and  is  a  function  of  only  the  electrostatic  (space- 
charge)  potential  of  the  plasma,  <f>.  The  slow  part  of  the  index  of  refraction  determined 
the  dispersion  relation,  Eq.  (13),  and  the  fast  part  determined  the  source  current  for  the 
harmonics.  Saturation  of  the  N*’*  harmonic  occurs  after  a  detuning  length,  Ld^  =  ctj,^  , 
pven  by  Eq.  (30).  The  harmonic  amplitude  is  maximum  when  the  laser-plasma  interaction 
length,  i,  is  an  integer  multiple  of  .  The  ratio  of  saturated  power  in  the  harmonic 
to  that  in  the  incident  laser  is  given  by  Eq.  (34).  This  expression  is  valid  for  long  lasers 
pulses,  ct£,  ^  \p  (no  10^®  cm~®  for  tx,  ~  1  psec),  and  for  interaction  distances  short 
compared  to  the  diffraction  length,  L  •C  Relativistic  harmonic  generation  favors  the 

use  of  high  densities  and  intense  lasers,  do  ~  1.  The  saturated  power  given  by  Eq.  (34), 
in  the  limit  dg  C  1,  is  a  factor  of  smaller,  where  6jv  are  constants,  than 

that  predicted  by  the  simplified  quiver  model,  which  neglects  the  self-consistent  plasma 
response.  This  reduction  in  the  harmonic  power  by  collective  plasma  effects  is  supported 
by  recent  experiments  on  harmonic  generation  in  pre-ionized  plasmas  [2],  [17],  [18]. 

The  effects  of  a  diffracting  incident  laser  field  with  a  finite  spot  size  ro  have  been 
analyzed  in  the  limit  dg  1.  Diffraction  is  important  for  interaction  lengths  L  >  itZr 
>  x-Zr  implies  S  =  2Ajprg/3  <  1  for  the  third  harmonic).  It  is  shown  that  no  third 
harmonic  signal  emerges  from  a  plasma  of  near  infinite  extent.  A  finite  third  harmonic 
signal  requires  the  use  of  a  semi-infinite  or  finite  slab  plasma.  The  third  harmonic  power 
emerging  from  the  edge  of  a  semi-infinite  plasma,  which  corresponds  to  the  focal  point  of 
the  incident  laser,  is  given  by  Eq.  (55)  and  is  a  factor  of  15  smaller  than  the  corresponding 
ID  saturation  power.  No  second  harmonics  are  generated  from  an  initially  uniform  plasma. 
The  generation  of  even  harmonics  requires  the  existence  of  transverse  gradients  in  the  initial 
plasma  density.  Circularly  polarized  light  will  also  generate  both  even  and  odd  harmonics 
when  initial  transverse  density  gradients  are  present. 
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The  most  severe  constraint  on  the  production  of  coherent  relativistic  harmon’c  radia¬ 
tion  is  that  of  phase  detuning.  In  a  dense  plasma  with  no  =  10*°  cm“®  (Ap  =  3.4  /xm)  and 
a  Ao  =  1  fim  laser  with  Jq  =  6.7  x  10^*  W/cm*  (ho  =  2.2),  the  saturation  efficiencies  for  the 
third  and  fifth  harmonics  are  P3IP0  =  2.2  x  10“®  and  Ps/Po  =  4.6  x  10"^°,  respectively. 
At  such  a  high  density  the  detuning  lengths  are  extremely  short,  =  8.0  /xm  for  the 
third  harmonic  and  Xj,  =  4.3  fim  for  the  fifth  harmonic.  If  a  scheme  for  phase  matching 
could  be  conceived,  the  interaction  distance,  L,  and,  hence,  the  harmonic  power,  Pn  ~  , 

could  be  dramatically  increased. 
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Fig.  1:  Schematic  of  an  intense  laser  interacting  with  a  plasma  to  produce  coherent  relativistic 
harmonic  radiation  in  the  forward  direction. 


27 


Fig.  2:  The  amplitude  of  the  normalized  vector  potential  of  the  harmonic  radiation,  jd,,  |,  for 
the  first  three  harmonics,  N  =  3,5,  and  7,  verses  the  normalized  propagation  distance, 
ct/Xo,  for  a  Ao  =  1  /im  wavelength  laser  of  intensity  Jo  =  6.7  x  10^*  W/cm*  (do  =  2.2) 
interacting  with  a  plasma  of  density  uq  =  10*®  cm  *  (Ap  =  3.4  /xm). 
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Fig.  3:  The  function  J7(j,ao),  proportional  to  the  third  harmonic  power  emerging  from  a 
plasma  of  width  2zo  centered  about  the  laser  focus,  verses  ao  =  2o/^Rt  for  6  = 
2jfepr2/3  =  10.0  (solid  curve),  1.0  (dashed  curve)  and  0.1  (dotted  curve),  where  Zr  = 
ibori  /2  is  the  Rayleigh  length. 


